On The Geometrical Description of Dynamical Stability II 
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Geometrization of dynamics using (non)-affine parametization of arc length with time is inves- 
tigated. The two archetypes of such parametrizations, the Eisenhart and the Jacobi metrics, are 
applied to a system of linear harmonic oscillators. Application of the Jacobi metric results in pos- 
itive values of geometrical lyapunov exponent. The non-physical instabilities are shown to be due 
to a non-affine parametrization. In addition the degree of instability is a monotonically increasing 
function of the fluctuations in the kinetic energy. We argue that the Jacobi metric gives equivalent 
results as Eisenhart metric for ergodic systems at equilibrium, where number of degrees of freedom 
N — > co. We conclude that, in addition to being computationally more expensive, geometrization 
using the Jacobi metric is meaningful only when the kinetic energy of the system is a positive 
constant. 



PACS numbers: 

Geometrization of dynamics is an alternative for sta- 
bility analysis of Hamiltonian systems, when normal lya- 
punov exponent is not[l|,|2|]. Its non-perturbative nature 
has opened the possibility for studying general condi- 
tions under which the traditional tools [3|, suitable for 
quasi-integrable systems, cannot be applied due to the 
strong mixing terms among the different degrees of free- 
dom. Moreover, it seems to be a promising framework 
for setting the statistical mechanics on dynamical and 
geometrical grounds. The general idea is to map the real 
motions q l (t) of the system 
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Pi = a ik (q)q k 



as geodesies on a properly defined manifold. Then, the 
quantitative stability analysis is obtained by means of the 
evolution with respect to arc length of the vector field 
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of the geodesic spread C( s ) 

is a congruence of geodesies, s is the arc length and r 
parametrizes the different geodesies. The general expres- 
sion for £ dynamics is given by the Jacobi-Levi-Civita 
(JLC) equationjj] 
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Where and Rj km are respectively the covariant deriva- 
tive along the geodesic and the Riemann-Christoffel ten- 
sor. Equation J2Jl, after opening the covariant deriva- 
tives, and taking into account the general expression for 
the geodesic equation 
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becomes 
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From the integration of Eq0J the geometrical indicator 
of stability is obtained: 



A = lim -log . , 
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Two widely used frameworks entail the application of 
the Eisenhart and the Jacobi metrics; the latter is given 
by (gj)ij = 2[E — V(q)]a,ij(q) and it is only applicable 
for time independent systems 0, 0. 

Usage of the Jacobi metric in the past has given the 
surprising suppression of chaos with increasing number 
of degrees of freedom 7] . This was linked to the approxi- 
mations done in the formalism; furthermore, parametric 
resonance was considered to be the fundamental source 
of instability arising from the fluctuating curvature. In a 
previous work |8( we addressed these controversial results 
aiming to warn about which conditions must be fulfilled 
in order to expect consistent results using either of the 
two metrics. In the context of a two dimensional sys- 
tem, we argued that the non-affine parametrization of 
arc length (s) with time (t) can be a significant source of 
non-physical instability. The danger of non-affine param- 
eterization seems to reside in the fact that trajectories 
are not compared at the same time and it is manifested 
geometrically by the large amplitude oscillations of the 
curvature. Moreover, from a mathematical point of view, 
it has been shown that due to the degeneracy of the Ja- 
cobi metric in the boundary set, there are geodesies that 
do not correspond to physical motions 0- 

In this letter we aim to make the consequences of 
the non-affine parametrization transparent and precise. 
Therefore, we use a well known paradigm of stability: 
a system of independent harmonic oscillators. We apply 
the stability analysis using Jacobi metric; any evidence of 
chaos in this dynamical system can then be understood 
as a failure of the method. Specifically, we investigate 
the dependence of the results obtained for the indicator 
of (in)stability on the fluctuation of the kinetic energy 
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along the geodesies. In contrast we will show analytically 
that the stability results in the framework of Eiscnhart 
metric are in accord with physical expectations. 
The time-independent Hamiltonian 



H(q,p) = 
V 



1 

2 

= S ik q 
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is our basic model. The solutions to the equations of 
motion are 



q k (t) = C k cos(ojt + i 



(7) 



where C k and 9 k depend on the initial conditions and 
in our case were chosen as C k — 1 and 9 k — k^jf, 
k = 1, ...,N. Here N stands for the number of degrees 
of freedom and the phases, 9 k , are homogeneously dis- 
tributed on a fraction (/) of 2ir (phase circle). Using Eq. 
O we get for the kinetic energy T and the square of its 
fluctuation 
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These settings allow us to tune the fluctuation y/a, vary- 
ing the values of / and N. It can be easily seen from the 
above equations that a decreases with increasing N, hav- 
ing the limit I S1 "i 27r {) I as N — > oo. When all of the oscilla- 
tors are in phase (/ = 0) , the variance takes its maximum 
value and the kinetic energy reaches the zero value every 
At = —. The specific dependence of a on the fraction 
of the phase circle for N = 10 is shown in FIG^ This 
is of interest under statistical considerations, where one 
expects the fluctuations of the kinetic energy to decrease 
as one makes more use of the phase circle. Generally in 
an TV-body bounded system, the kinetic energy is a fluc- 
tuating quantity. In non-integrable systems due to the 
mixing terms (non-separability), these fluctuations are 
normally reduced with increasing number of degrees of 
freedom. Since our toy model has no interactions at all 
and is coherent, it is additive but not necessarily ergodic; 
therefore a does not vanish when N — > oo. However, our 
main interest is to explicitly show the shortcomings of 
the Jacobi metric, which is unaffected by the properties 
of our toy model. Moreover, we will show that A strictly 
depends on a. 

Eq[3 when applied to EqJBJ using Jacobi metric, be- 
comes 
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This is radically different (due to the extra mixing terms) 
to that obtained from Eiscnhart's metric, which corre- 
sponds to the tangent dynamic equations of motion 
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with the couplings /■ , J k and K k given by 
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It is perhaps worth while mentioning that the solution 
of Ea. l|ll|l evidently gives A < for any setting of initial 
conditions; moreover A = is of measure zero. Thus, 
the Eisenhart geometrization behaves properly for this 
system and is a good indicator of (in)stability. Once the 
analytic solutions Q were substituted in Eq. fTTijl it takes 
the final form 



d 2 e 2, k 
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dt 



[J k + K k ]e=0 (12) 



The above expressions have a typical time dependence 
and structure from which instabilities such as parametric 
resonance may arise with the extra ingredient that 
the amplitude of the couplings are generally unbounded 
due to the inverse powers of T. Although these expres- 
sions are system dependent, the appearance of inverse 
powers of T is a general feature of the Jacobi metric and 
it is responsible for the singularities in the boundary set. 

In order to make our analysis quantitative, A was 
obtained from the numerical integration of £ dynamics 
Eci. (|12|) computed using different values of / and N. In 
FIG|21 we superimpose these values of A as functions of 
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FIG. 1: The polar plots for the square of the fluctuations in 
the kinetic energy a vs. the fraction of the phase circle / for 
N = 10, u = 2tt. 



yfa (i.e. fluctuations in the kinetic energy), obtained 
from the EqOH for different settings. The fact that the 
plot shows a smooth dependence implies that the rel- 
evant quantity of interest, on which A depends, is the 
fluctuation of the kinetic energy. 

Plot of X vs. The fluctuations in T 
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FIG. 2: The superposition of the geometrical lyapunov ex- 
ponent A vs. the fluctuations of the kinetic energy i/a for 
TV = 2 + f and / = 0.05,0.1, ...,0.45, where j = 1,..., 14. 

In the light of a previous work8j, FIG |3 makes ev- 
ident the strong correlation between the values of the 
geometrical stability indicator A and the variance in the 
kinetic energy with initial conditions for the two oscil- 
lators. Needless to say, strong non-physical divergences 
arise from the "kicks" produced in the ^-dynamics when 
the kinetic energy approaches zero periodically. As we 
discussed before, these "kicks" are nothing but diver- 
gences of the curvature tensor close to the boundary set 
because the curvature tensor contains inverse powers of 
the kinetic energy. Lastly we show, for the maximum 
uniform distribution, i.e./ = 1, the dependence of A on 



FIG. 3: Comparison of A and the absolute variance for two 
degrees of freedom vs. /. 



fundamental frequency of the oscillators lu. As is evi- 
dent in the FIG 01 the lyapunov exponent increases with 
increasing uj. Curiously, although a — according to 
Eq|jJ we obtained small but positive A. The reason is 
that equations like (|12f) have extremely unstable behav- 
ior. Therefore, even in cases where the variance was set 
to zero, small fluctuations produced by numerical evalua- 
tion are amplified, leading to positive exponents, that are 
naturally reduced with reducing the value of (u>), which 
multiplies the mixing terms in Ea ll2l 
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FIG. 4: The geometrical lyapunov exponent A vs. the fun- 
damental frequency of oscillation u) for N = 10 and / = 1. 

In this letter we have shown that in the framework of 
the Jacobi metric the evolution of the vector field of the 
geodesic spread is not physically meaningful, ^-dynamics 
shows an exponential divergence A > 0, for a system 
which is physically stable. We showed that these non- 
physical divergences are due to the fluctuations in the 
kinetic energy, which for ergodic systems should vanish at 
the thermodynamical limit. In addition these oscillations 
have the cumulative effect of parametric resonance-type 
that can be confused with real instability. This instabil- 
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ity was induced by the non-affine parametrization of the 
arc length with time which causes the local comparison 
of the trajectories to be at different times. It is expected 
that, only for ergodic systems at the limits of N — > oo and 
t — ► oo j the Jacobi metric would give equivalent results 
as the Eisenhart metric would. Moreover, the stability of 
this system of independent oscillators is shown to be de- 
pendent on the fundamental frequency which of course 
is incorrect. Instabilities show up even in cases where 
the variance of T is theoretically zero by the choice of 
the initial conditions. It is an evidence of the unstable 
character of the solutions of the JLC equation for the 
Jacobi metric, which amplify the error introduced by nu- 
merical integration. Therefore, the Jacobi metric gives 
correct stability results only when the kinetic energy of 
the system as whole is strictly a positive constant, which 
is practically impossible. 
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